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"Abstract 

[We study the dynamics of the many-body state of ultracold bosons trapped in a bistable optical lattice in an optome- 
chanical resonator controlled by a time-dependent input field. We focus on the dynamics of the many-body system 
[following discontinuous jumps of the intracavity field. We identify experimentally realizable parameters for the bistable 
quantum phase transition between Mott insulator and superfluid. 
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> '1. Introduction 

^ , The coupling of coherent optical systems to micromc- 
^H'chanical devices, combined with breakthroughs in nanofab- 
I— ixication and in ultracold science, has opened up an cxcit- 
ing new field of research, cavity optomechanics. Several 
^ groups have now demonstrated very significant cooling of 
the vibrational motion of a broad range of moving mir- 
■rj" rors, from nanoscalc cantilevers to LIGO-class mirrors [ij, 
cind there is every reason to believe that the quantum me- 
chanical ground state of motion of these systems will soon 
be achieved. In a parallel development, ultracold gases 
as well as Bose-Einstein condensates placed inside optical 
•resonators have been shown to behave under appropriate 
conditions much like moving mirrors Following these 

developments, cavity optomechanics is rapidly becoming a 
yery active sub-field of fundamental and applied research 
cit the boundary between AMO physics, condensed matter 
5^ physics, and nanoscience. 

Cavity optomechanics presents considerable promise 
both in opening the way to address fundamental ques- 
tions related to pushing quantum mechanics toward in- 
creasingly macroscopic systems, but also in applications 
that span a variety of areas from quantum detection to 
the coherent control of microscopic atomic and molecu- 
lar systems and/or of nanoscale devices. Specific exam- 
ples include nanomechanical cantilevers coupled to a Bose- 
Einstein condensate [3| or to dipolar molecules , to a sin- 
gle atom 0, oscillating mirrors coupled to atomic vapors 
8|, magnetized resonator tip to color center of diamond 
9(] . Cold atoms can probe the state of the nanostructure, 
or conversely, the optomechanical cavity setup can serve as 
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a diagnostic tool for the quantum state of atoms trapped 



inside the resonator [ll|, . On a more applied side such 



systems will enable the development of ultrasensitive force 
sensors and may find applications in quantum information 
processing technology. 

In recent work we studied the many-body state of ul- 
tracold bosons in a bistable optical lattice potential in 
an optomechanical resonator. We considered explicitly 
the weak-coupling limit where the coupling between the 
cavity-field and the movable mirror results in a bistable 
optical lattice potential for the atoms, and showed how 
such a cavity plus cold-atom system can be engineered so 
that a superfluid and a Mott insulator phase are bistable 
ground states for the ultracold atomic gas [l3l |. 

The present paper extends this study to consider the 
dynamics of the transition between these two states as 
the optical potential undergoes a bistable loop. A time- 
dependent variational principle combined with a Gutzwiller 
ansatz arc introduced to calculate the evolution of the 
quantum many- body state as a super-Gaussian incident 
light pulse [ij] switches the intracavity optical field. We 
discuss conditions under which this process can proceed 
adiabatically, and identify experimental parameters that 
would permit to test our predictions. 

Section II reviews our model of a quantum-degenerate 
Bose gas confined by the standing wave generated in a 
Fabry-Pcrot cavity with one movable end mirror, and briefly 
reviews its steady-state properties. Section III introduces 
an effective potential description of the system and com- 
ments on the conditions required for an adiabatic evolution 
of the atoms. It then turns to the time evolution of the 
superfluid order parameter during the bistable superfluid 
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driven by radiation pressure, 



Figure 1: A Fabry-Perot cavity with a movable mirror displaced 
from its equilibrium position by q. 
cavity and transmitted light amplitudes, respectively. 



to Mott insulator transition. Finally, Section IV is a con- 
clusion and outlook. 



2. Model 

We consider a Fabry-Pcrot cavity with one fixed end 
mirror and the other one, of mass M , mounted on a spring 
of frequency and damping rate 7™. An optical field 
Ep of frequency ujp is incident on the fixed mirror. The 
reflectivity of both mirrors is R, with i? ~ 1 so that the 
intracavity field can be approximated by a standing wave 
of amplitude Ec- In the absence of that field the movable 
mirror is at its equilibrium position q = and the cavity 
length is Lq, see Fig. [TJ 

A sample of ultracold bosonic two-level atoms with 
transition frequency LOa is loaded inside the optical lattice 
formed by the standing wave, and the atoms are assumed 
to interact with the light field in the weak coupling regime 
NgQ/\Aa\ <^ Kc, where N is the total number of atoms, 
go is the single-photon atom-field coupling strength, Kc is 
the cavity decay rate, and = lUp — uja is the atom- 
field detuning [15[. In this limit the intracavity field has 
no significant dependence on the atomic distribution and 
we can investigate the dynamics of the intracavity field by 
using the theory of an empty cavity. In this paper we de- 
scribe both the light field and the movable mirror as clas- 
sical objects, while the center-of-mass motion of the atoms 
is quantum mechanical. The extension to mirrors cooled 
near their quantum mechanical ground state of vibration 
will be considered in future work. 

Ignoring the retardation effects due to the propagation 
of the light field back and forth inside the cavity, the intra- 
cavity field is governed by the familiar equation of motion 

dEc , . . , „ 2k, 



dt 
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where Ac = u)p — ujc is the cavity-pump detuning and T 
is the transmittivity of the mirrors, with i? -I- T = 1 in 
the absence of mirror absorption. The movable mirror is 
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where Frp ~ AeQ\Ec\'^ / A is the radiation pressure force, 
with A being the cross-sectional area of the laser beam. 

Equations ((T|) and ^ are coupled through the cavity- 
pump detuning since the cavity frequency is g-dependent. 



q 



(1 - q/Lo) 



(3) 



where n is the integer associated with the mode closest to 
the laser frequency in the single-mode theory considered 
here. The approximate equality assumes that the mir- 
ror displacement due to radiation pressure is small com- 
pared to both Lq and to the cavity mode wavelength Ac = 
2Trc/uJc- 

It is convenient at this point to introduce dimensionless 
units and to scale times to fl"^, lengths to Ap/2 and elec- 
tric field amplitudes to y^2Er/ca€Q, where E^ is the pho- 
ton recoil energy and a — 37rc^7a/(2a;j^Ao), with the 
natural linewidth of the atomic resonance. Equations ([T]) 
and ^ reduce then to the dimensionless form 
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where r is the (dimensionless) time, X and Y the in- 
tracavity and incident field amplitudes, ^ mirror's dis- 
placement, K and 7 the intracavity optical field and mov- 
ing mirror damping rates, and E the offset from the cav- 
ity resonance in the absence of radiation. Finally, /3 = 
{AEr)l[\pMVi?-ca). 

The quantum-degenerate atomic sample, assumed to 
be at zero temperature is trapped in the optical lattice 
along cavity axis, its transverse motion being strongly con- 
fined by a strong transverse potential. Its Hamiltonian can 
be then approximated as the one-dimensional form 



dx^'^{x) 



where 4' (a;) is the Schrodinger field operator, g is the two- 
body interaction coefficient, m is the atomic mass and 
Vb = OLCtQ\Ec^ jl = \X^Er is the depth of the optical 
potential. In the tight-binding approximation this many- 
body Hamiltonian can be simplified to a single-band Bose- 
Hubbard Hamiltonian 



^BH = -J^ajaj 



1) 



(6) 



where hi is the bosonic annihilation operator for site z, n, 
is the corresponding number operator, the subscript 
labeling nearest neighbor pairs. Finally J is the inter-site 



2 



•120 
a 
B 
a 



> 

03 

o 



T3 
in 

O 
C/2 



10 









1 


IF 





0.03 0.04 0.05 

Scaled Input Intensity |Yp 



0.06 



Figure 2: Bistable dimcnsionless intracavity intensity p as a func- 
tion of the dimcnsionless pump intensity |yp for T = 0.01, initial 
offset 5 = —0.0035 and (3 = 0.0002. In real units, and for the pa- 
rameters considered here, the threshold intensities for the bistable 
domain are 0.74mW and 0.87mW. The steady-state position of the 
mirror is likewise bistable. 



tunneling matrix element and U is the pair interaction 
energy. The parameters J and U can be evaluated by ex- 
panding the field operators on the Wannicr basis of the 
lowest Bloch band, and then evaluating the pertinent in- 
tegrals [l^ . 

Setting the time-derivatives to zero yields the steady- 
state solution 



2Y/VT 



i-i2TT + 
P\XA'- 



(7) 



Eliminating the equilibrium displacement f s in these equa- 
tions results in a transcendental equation for the intracav- 
ity field intensity iX^p, 



4/r 



(8) 



This equation is known to have a bistable solution for an 
appropriate choice of parameters, see Fig. [51 The upper 
and lower branch with positive slope are generally stable, 
while the dashed branch with negative slope is unstable 
[l9| . It is clear from Eq. ([7|) that the steady-state dis- 
placement of the movable mirror has the same bistable 
property as the intracavity field intensity. 

As is well known, the ground state of the atomic system 
is determined by the ratio J/U which can be controlled by 
varying the lattice depth Vq. For a shallow potential, in- 
terwell tunneling dominates and the many-body ground 
state is a supcrfluid, while for a deep enough potential on- 
site interactions dominate and the atoms enter the Mott- 
insulator phase with integral filling factor [l3|. The dot- 
ted line in Fig. [2] indicates the critical depth Vq at which 



'^AU calculations presented in this paper are for -^^Na and for 
Ap=985nm, M= 0.078g, Q, = 2it x lOHz and mirror reflectivities 
R = 0.99. 



the superfiuid-Mott insulator phase transition occurs, for 
sodium atoms trapped by a laser of wavelength 985nm. Of 
particular interest is the bistable region where the shallow 
lower branch corresponds to a superfiuid phase and the 
deeper upper branch to a Mott insulator phase. The next 
section discusses the dynamics of the transition between 
these two phases as the optical field is varied along the 
bistable loop. 

3. Dynamics 

Three important time scales are relevant in understand- 
ing the dynamics of the system: the characteristic time 
of the mirror, the cavity buildup time and the 
interwell tunneling time h/J of the atoms. In the weak 
coupling limit the dynamics of the cavity system (the light 
field plus the movable mirror) is determined by the former 
two. Furthermore, in the case of a bad cavity the dynamics 
of the light field is much faster than that of the movable 
mirror, ^ so that it can adiabatically follow the 
instantaneous displacement of the movable mirror. Note 
that this approximation neglects the non-adiabatic effects 
that lead to optical damping, see Ref. [2^ and references 
therein. 

This paper treats the motion of the moving mirror clas- 
sically as it is not concerned with the dynamics of mirror 
cooling. We therefore ignore optical damping in the fol- 
lowing, and adiabatically eliminate the optical field in the 
Eq. ([5]), resulting in the nonlinear oscillator equation 



4\Y\yT 
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To further analyze the dynamics of this nonlinear os- 
cillator we introduce the effective potential 21 1 



?0 



4|y|2/r 



co 1 + 47r2 (,5 + ey /r^ 



de (10) 



where the first term corresponds to the restoring force on 
the movable mirror and the second to the radiation pres- 
sure force. Figure [3] shows that potential for three values 
of the input intensity. For a weak input field, Fig.jJlJa) the 
potential has a single minimum, but increasing it sees the 
appearance of a second local minimum. Fig. [31(b), indica- 
tive of bistability. Further increasing the input intensity 
past the bistable region the initial minimum disappears as 
expected, see Fig. |3Ic). 

For concrctencss wc assume in the following that the 
optical field incident on the Fabry-Perot is a super-Gaussian 
pulse. Fig. mja). Since the intracavity field is assumed to 
follow the motion of the mirror adiabatically, this readily 
yields the time-dependent optical potential 



|X(r)|^ 



4|y(r)|2/r 
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Figure SKb) illustrates the dynamics of the intracav- 
ity intensity as the incident intensity is varied to switch 
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Figure 3: Potential V{i) for (a) input intensity |Yp = 0.03; (b) 
|yp = 0.044; (c) |yp = 0.06. Other parameters are the same as in 
Fig. 2. All variables are dimensionless. 
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Figure 4: (a) Input intensity |yp as a function of dimensionless 
time T. The two dashed lines label the critical intensities at the 
beginning and the end of the bistable region, (b) Evolution of the 
intracavity intensity |Xp controlled by Here 7 = 0.5 and the 

other parameters are as in Fig. 2. The dotted line labels the critical 
depth of the atomic phase transition. 



the system from the lower to the upper bistability branch 
and back. In this example |y(r)p is initially switched to 
a value only slightly above the critical value \Ya^ where 
the first local minimum of the effective potential degener- 
ates into a plateau, hence it takes a relatively long time to 
switch to its new steady-state value. Under appropriate 
circumstances this critical slowing down (2^ can be ex- 
ploited to allow the atoms to adiabatically follow the po- 
tential changes in the superfluid region, see below. Follow- 
ing that slow evolution stage, |X(r)p grows nearly expo- 
nentially as the mirror falls into the newly formed potential 
well and completes its transition to the upper branch. For 
the switch back to the lower branch shown in this exam- 
ple, the lower value of |y(r)p is chosen to be significantly 
below lYbp, so as to reduce the effects of critical slowing 
down. As expected, |X(r)p then switches much faster to 
the lower branch. 

The time dependence of the potential |X(t)P results in 
Wannicr functions that are likewise time-dependent, and 
hence also in time-dependent tunneling rate J(r) and on- 
site interaction U{t). In order for the many-body ground 
state of the atoms to follow adiabatically the changes in 
potential, two requirements must be fulfilled. The first one 
is that the variation of the potential depth must be slow 
enough to prevent the occurrence of inter-band excitations, 
that is the atoms must remain in the first Bloch band at all 
times. This condition is usually easy to satisfy for atoms 
with quasi-momentum g ~ because of the existence of a 
band gap [23|, |24| . The explicit adiabaticity criterion is 



at 

We note that for the experimental parameters considered 
here this condition is satisfied for even the fastest tran- 
sients in Fig. 4(b): The maximum value reached by \ ^Vo/Er\ 
is 2 r2 which is much less than 16 ujr- We therefore assume 
in the following that that inter-band excitations are negli- 
gible, and evaluate J{t) and U (t) with Wannicr functions 
of the first band. 

In addition to this single-particle adiabaticity condi- 
tion, we also need to consider the time scale associated 
with many-body effects. As the potential depth varies in 
time, the atoms need a sufficient amount of time to tunnel 
and redistribute themselves across the lattice, and hence 
to settle in their new ground state — think specifically of a 
transition between a Mott insulator and a superfiuid. Al- 
though this is not a significant consideration if the atoms 
remain in the Mott insulator phase, which is quite insen- 
sitive to inter- well rearrangement, it is more important in 
the superfluid region, due to the variation of the tunnel- 
ing rate J(t)[2^, [2^. That is the reason why the critical 
slowing down can be useful in controlling the variation of 
the potential in the superfluid region. Whether or not 
the second adiabaticity condition is satisfied can be deter- 
mined from the time-dependent order parameter which we 
calculate next. 
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The evolution of the atomic ground state can be ob- 
tained from the time-dependent variational principle 27, 

S{G\th^^-HBB{t)\G) ^0, (12) 

where the Bose-Hubbard Hamiltonian ifsH (t) depends on 
time via the coefficients J{t) and U{t), and \G) is taken 
to be given by the time-dependent mean-field Gutzwiller 
ansatz 



\G) 




(13) 



Here Ni is the number of lattice sites 

c{i) 



are Fock states, 

and fn'{t) are probability amplitudes that preserves the 
normalization of the wave function. We assume that the 
lattice is uniform with unit filling factor and insert the 
Gutzwiller ansatz into Eq. p^ . This yields the set of 
coupled nonlinear equations for the amplitudes /n(t) 



-2J{t) {V^fn-i (a) + (at) 
Uit) 



-n {n + 1) /„ 



(14) 



where (a) = X^J^o Vn+Tf*fn+i is the atomic superfluid 
order parameter. Our numerical results are for an initial 
optical potential depth of 5-B,-, for which the tight binding 
approximation is valid and the many-body ground state 
is superfluid. Figure [D^a) shows the ratio 2J{t)/U{t) as 
a function of time. The transition between the super- 
fluid and Mott insulator phases is clearly shown by the 
absolute value of the superfluid order parameter | (a) | plot- 
ted in Fig. [5]^b). As the ratio 2J/U decreases below the 
critical value 0.17, corresponding to the critical potential 
depth Vq = lO.SEr for sodium atoms, the order parame- 
ter rapidly drops to zero, indicating that the atoms enter 
the Mott insulator phase. When 2J/U rises back above 
that critical value, the order parameter recovers a non- 
zero value, indicative of the return of the atoms to the 
superfluid phase. 

The small oscillations of the order parameter in the 
Mott-insulator region are due to collapses and revivals of 
the condensate with period T = h/U [291, i30|. The oscilla- 
tions following the return of the system to the superfluid 
region were first predicted in Ref. [sij. They result from 
a rapid transition of the system from being a Mott insu- 
lator to a superfluid, and can be suppressed via a slower 
variation of potential depth. The small amplitude of these 
oscillations, less than 5 percent of the initial value of the 
order parameter, indicates that the second adiabatic con- 
dition (associated with tunneling) is well satisfied in this 
example. 

For comparison Fig. 5(c) shows the time-dependent su- 
perfluid order parameter for a situation here non-adiabatic 
effects are more apparent. In this case the mass and vibra- 
tion frequency of the mirror have been changed to O.OSlg 
and 2tt x 50Hz. Since these parameters give the same 
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Figure 5: (a) 2J{t)/U{t) as a function of the dimensionless time 
T. The dotted line indicates the critical value 2J/U = 0.17, corre- 
sponding to the potential depth lO.SEr for our parameters, for the 
superfluid and to Mott insulator phase transition, (b) Dynamics 
of the superfluid order parameter |(a)| during the transition from 
the superfluid to the Mott insulator and back to the superfluid. 
Here \p = 985 nm and uJr = Er/h = 2it X 8.9 kHz for ^^Na. The 
mass and vibration frequency of the mirror are M = 0.078 g and 
f2 = 27r X 10 Hz, respectively. In real units, the full cycle through the 
bistable loop takes 2.5s. (c) Same as (b) but for M = O.OSlg and 
Q = 27r X 50Hz 
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steady-state solutions as Fig. 2, the form of the evolution 
of the intracavity intensity is unchanged, but it occurs now 
on a shorter time scale. The cycle through the bistable 
loop takes now 0.5s in real units, five times faster than the 
case of Fig. 5(b). Although the first adiabatic condition 
is still satisfied for these parameters, it is apparent from 
Fig. 5(c) that the amplitude of the oscillations in both the 
Mott insulator and the superfluid regions is much larger 
than in Fig. 5(b) an indication of the breakdown of the 
second adiabatic condition. 

We conclude by remarking that the finite temperature 
of the classical mirror leads to thermal position fluctua- 
tions which result in turn in fluctuations for the intra- 
cavity field strength the lattice potential depth SVq. In 
order to prevent this effect we request that SVq <^ Vnun, 
where Vmin is the minimum value of the potential depth 
during the switching process. This inequality gives a mir- 
ror temperature limit of T <C lO^K. 

4. Conclusion 

We have analyzed the dynamics of a quantum-degenerate 
ultracold sample of bosonic atoms trapped in the optical 
lattice of a bistable optomechanical system. The varia- 
tion of intracavity optical lattice depth results in a bistable 
quantum phase transition between a superfluid and a Mott- 
insulator ground state. We considered the concrete case of 
a super-Gaussian incident light pulse to drive the atomic 
system around the bistable loop, and found that the crit- 
ical slowing down of the optomechanical system can be 
useful to prevent the excitation of the atoms in a po- 
tential ramping-up stage. Our numerical results indicate 
that the bistable superfluid to Mott-insulator phase tran- 
sition should be realizable for realistic experimental pa- 
rameters. Future work will extend these studies to the 
strong-coupling regime, where the optical field depends on 
the motion of the cold atoms as well as the mirror. We 
will also investigate to which extent the mirror can serve 
as a measurement device to observe the dynamics of a 
quantum phase transition nondestructively, in particular 
in cases where it is cooled to near its quantum mechani- 
cal ground state of vibration. Additional goals include the 
study of the optically induced quantum correlations and 
entanglement between the ultracold atomic system and the 
mirror, the quantum control of the mirror motion by the 
atoms, and conversely of the atoms by the nanomechani- 
cal system, and finally the development of novel types of 
quantum sensors. 
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